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Abstract 



Damping in a dilute Bose gas is investigated using a mean-field approxima- 
tion which describes the coupled oscillations of condensate and non- condensate 
atoms in the collisionless regime. Explicit results for both Landau and Beli- 
aev damping rates are given for non-uniform gases. In the case of uniform 
systems we obtain results for the damping of phonons both at zero and fi- 
nite temperature. The isothermal compressibility of a uniform gas is also 

discussed. 
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I. INTRODUCTION 



The low-lying collective excitations of magnetically trapped Bose gases have been the 
object of very accurate experimental measurements The first experiments were 

carried out at low temperature, with approximately all the atoms in the condensate state. 
These experiments showed almost undamped oscillations of the condensate at frequencies 
which have been found in excellent agreement with theoretical predictions based on the 
T = time-dependent Gross-Pitaevskii equation In more recent experiments the 
study of the low-energy collective modes has been extended to higher temperatures where 
the condensate oscillates in the presence of a considerably large fraction of above- condensate 
atoms. In this case, evidence is given of large frequency shifts with unexpected features and 
of strong damping rates which have not yet been understood theoretically. 

On the theoretical side, extensions to finite temperature of the Gross-Pitaevskii equation 
have been put forward fSj-Q and have been very successful in explaining experimental results 
on the thermodynamic properties of these systems, such as condensate fraction, internal en- 
ergy, specific heat and critical temperature These mean- field descriptions, however, 
seem to be inadequate when applied to the study of collective excitations at finite tempera- 
ture [§,3. First of all, the proposed mean-field approximations do not account for damping. 
Secondly, the resonance frequencies are predicted to vary with temperature mainly because 
the number of atoms in the condensate changes, and the dependence of the frequency shift 
upon this quantity is expected to show the same behaviour as the corresponding T = 
dependence on the total number of atoms in the trap. None of the features exhibited by 
the experimental data on the frequency shift can be explained using these descriptions. 

In the presently available finite-T extensions of the Gross-Pitaevskii equation it is as- 
sumed that the condensate oscillates in a bath of thermally excited atoms at rest and in 
thermal equilibrium in the effective mean-field potential generated by the average conden- 
sate density. This assumption is valid if the time scale on which the thermal cloud oscillates 
is much larger than the inverse frequency of the condensate oscillations. Since both the 
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condensate and the thermal cloud vary on comparable time scales of the order of the trap 
frequency, a full dynamic treatment of the coupled oscillations of the two clouds should be 
developed. The dynamic effects we are aiming to describe are irrelevant in the calculation 
of the thermodynamic properties of the system, for which the standard finite-T mean-field 
approaches are well suited. In fact, as pointed out in Refs. [0,|Tni) the fine details of the 
excitation energies do not affect the thermodynamic behaviour, for which what matters is 
the density of states at a given energy. 

Another problem arises as to whether the appropriate regime to describe the experimental 
situation in [Q] is collisional or collisionless. In the coUisional regime, which takes place at 
high temperatures and densities, the damping mechanisms are of dissipative type and the 
dynamic of the system is described by two-fluid hydrodynamics, recently developed for 
trapped gases in Ref. |Tl|. For very dilute systems at low temperatures the mean free 



path of the elementary excitations becomes comparable with the size of the system and 
collisions play a minor role. Damping in this regime is not related to thermalization processes 
but to coupling between excitations, and can be described in the framework of mean field 
theories. As suggested in Ref. |T2|, the collisionless regime may be appropriate for the JILA 



experiments [^0, but probably not for the MIT experiments 0. 

In the collisionless regime and at finite temperature the damping of the low-lying col- 
lective modes is dominated by Landau damping. The idea that this mechanism might be 
relevant to explain the damping rates in trapped Bose gases was first suggested by Liu 
and Schieve |jl3[ and has then been developed in Refs. |1T^-|16|. In Ref. Pitaevskii and 



Stringari have derived, using perturbation theory, an expression for Landau damping which 
is applicable to trapped Bose gases, and have shown that, applied to the uniform case, it 
reproduces known results for both the low- and high-temperature asymptotic behaviour of 
the phonon damping. 

In this paper we develop a time-dependent mean-field approach based on the Popov 
approximation, which describes the dynamics of a Bose-condensed gas in the collisionless 
regime. We obtain a set of coupled equations for the condensate and non-condensate com- 



ponents which allow us to calculate damping coefficients. In non-uniform gases we derive 
explicit expressions both for Landau damping, which coincides with the finding of Ref. ||14|| , 



and for Beliaev damping. In the uniform case we reproduce all known results on the damping 
rates of phonons both at T = and at finite T. 

However, the calculation of the temperature dependence of the frequency shifts using 
the Popov approximation is not reliable. For a uniform system the calculation should give, 
in the long wavelength limit, the correction to the velocity of zeroth sound due to quantum 
and thermal fluctuations. On the other hand, at T = 0, the velocity of zeroth sound 
is directly related to the bulk compressibility of the system, for which we show that the 
Popov approximation gives an incorrect result. By studying the isothermal compressibility 
we also find that the Popov approximation is inconsistent at low temperatures, because it 
neglects fluctuations which are relevant for temperatures smaller than the chemical potential. 
Since the same fluctuations are also important in the calculation of the velocity of zeroth 
sound in the low temperature regime, we draw the conclusion that a dynamic mean-field 
description, which gives correct account of both damping rates and frequency shifts, can 
only be developed going beyond the Popov approximation. 

The paper is organized as follows: in section 2 we develop the formalism of the time- 
dependent mean-field approximation. In section 3 we study the damping of the oscillations 
in a non-uniform system obtaining explicit expressions for the Landau and Beliaev damping. 
In sections 4-A, 4-B we apply the results of section 3 to the uniform case. In section 4-C we 
investigate the isothermal compressibility of a uniform gas beyond the Popov approximation. 

II. TIME-DEPENDENT MEAN-FIELD APPROXIMATION 

In the presence of a non-uniform external field V^xtlr) the grand-canonical Hamiltonian 
of the system has the form 
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+ |y cir^t(r,t)^t(r,t)^(r,t)^(r,t) (1) 

in terms of the creation and annihilation particle field operators '?/'^(r, t) and ip{r,t). In the 
above equation g is the interaction coupling constant, which to lowest order in the s-wave 
scattering length a is given by g = Aixfi^a/m. The equation of motion for the particle field 
operator then follows immediately and reads 

ifi^^ipiv.t) = [V'(r,t),K] 

+ V;,i(r)-/i U(r,t) + g^t(r,t)^(r,t)^(r,t) . (2) 



Y 2m 

According to the usual treatment for Bose systems with broken gauge symmetry we define 
a time-dependent condensate wavefunction $(r,t) fl^ 

$(r,t) = (^(r,t)) , (3) 

which allows us to describe situations where the system is displaced from equihbrium and 
the condensate is oscillating in time. The average (...) in equation (||) is thus intended to be 
a non-equilibrium average, while time-independent equilibrium averages will be indicated in 
this paper with the symbol (...)o- The particle field operator can then be decomposed into 
a condensate and a non-condensate component 

^(r,t) = $(r,t)+V;(r,t) , (4) 

and the non-condensate term satisfies, by definition, the condition {ip{v, t)) = 0. By applying 
the decomposition (^) to the Heisenberg equation (Q) the term cubic in the field operators 
becomes (all quantities depend on r and t) 

-I- + . (5) 

We assume that for dilute systems the cubic product of the non-condensate operators (last 
term in eq. (^) has a negligible effect on the dynamics of the condensate and we set its 
average value equal to zero: 



(^t(r,t)^(r,t)^(r,t)) = . (6) 

One thus obtains the following equation for the time rate of change of {tp{r,t)) 

ih-<^{r,t) = ( + Ve.t{r)-fij $(r,t) +g|$(r,t)p$(r,t) 

+ 2g$(r,t)n(r,t)+g$*(r,t)m(r,t) , (7) 

where we have introduced the normal and anomalous time-dependent densities defined, 
respectively, as 

n(r,t) = (^t(r,t)^(r,t)) 

m{r,t) = {4j{r,t)^{r,t)) . (8) 

The equation of motion (|^) with n = rh = corresponds to the usual T = Gross-Pitaevskii 
equation for the condensate wavefunction, while its extension including the normal and 
anomalous densities has been already discussed by many authors both in the study of ther- 
modynamic properties and of the collective modes at finite temperature [^|-|9|,pJ]] . The nov- 
elty here is that we will treat both terms within a time- dependent mean-field approximation 
holding in the coUisionless regime which is the object of the present work [Q. 

We are interested in the small-amplitude regime in which the condensate is only slightly 
displaced from its stationary value $o(r) = {ip{r))Q 

$(r,t) = $o(r) + 5$(r,t) , (9) 

where (5$(r, t) is a small fluctuation. In the same way we consider small fluctuations of the 
normal and anomalous densities 

h{r, t) = n°(r) + 6h{r, t) 

m{r,t) = Th°{r) + 5m{r,t) (10) 

around their equilibrium values fi°(r) = {ip''{r)ip{r))o and m^{r) = {ip{r)ip{r))Q. 

In the so-called Popov approximation the effects arising from the equilibrium value of 
the anomalous density in eq. (|^) are neglected and the following ansatz is introduced in the 
mean-field scheme 
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m°(r) = . (11) 

This approximation has been first introduced by Popov in the study of a uniform weakly 
interacting Bose gas at finite temperature (for a detailed discussion see Ref. and 



). More recently the Popov approximation has been extensively used in the study of 
properties of magnetically trapped Bose gases at finite temperature [PH7|.|9|,p| . The Popov 
approximation gives a gapless spectrum of elementary excitations which at T = coincides 
with the well-known Bogoliubov dispersion relation, while at high T it approaches the finite- 
temperature Hartree-Fock spectrum ||2T[| . 

By using the Popov prescription (|TTD, the real wavefunction $o(i") satisfies the stationary 
equation 

+ VUv) - /i + g(no(r) + 2n0(r))^ %{v) = , (12) 

where r;,o(r) = |$o(r)P is the condensate density, while the time-dependent equation for 
5$(r, t) is obtained by linearizing the equation of motion (|^) 

ihj^5^{r, t) = + Ve.t{r) - /i + 2gn(r)^ <5$(r, t) + gno(r)5$*(r, t) 

-f2g$o(r)5n(r,t)+g<l>o(r)5m(r,t) , (13) 

where we have introduced the total equilibrium density n(r) = no(r) + n°(r). 

From equation (0) one clearly sees that the oscillations of the condensate are coupled 
to the fluctuations 5fi(r, t) and 5m(r, t) of the normal and anomalous densities for which we 
need the equations of motion. To this aim it is convenient to express the non-condensate 
operators -ip, ip'^ in terms of quasiparticle operators aj, a] by means of the Bogoliubov linear 
transformations 



V^t(r,t) = Y: {u*ir)a}{t) + v,{r)a,{t)) . (14) 
j 

The normalization condition for the functions Uj{r), Vj{r), which ensures that the quasipar- 
ticle operators aj, a] satisfy Bose commutation relations, reads 



J dr {u*{r)uj{r) -v*{r)vj{r)) = 5ij 



(15) 



By using the transformations (|l^) the quantities n{r,t) and m{r,t) can be expressed in 
terms of the normal and anomalous quasiparticle distribution functions defined by 



hit) = {al{t)aM 
Qijit) = {ai{t)aj{t)) . 

The time evolution of these functions is fixed by the following equations of motion 



(16) 



a I 

d , , 



a\{t)a,{t),K 
ai{t)aj{t), K 



(17) 



In order to calculate the commutators of eq. (0) one notices that, after substituting the 
decomposition (^ of the particle field operator into the Hamiltonian (P, only the terms 
quadratic and quartic in the non-condensate operators ip, ip'^ give non- vanishing contribu- 
tions, because, according to eq. (|]), we set to zero all averages of cubic products of the 
non-condensate operators. The terms in the grand-canonical Hamiltonian relevant for the 
calculation of the commutators are thus given by 



K2 + Ki 



^^(r,t)| 



2m 



+ Ve,t{r) - /i 1 ^(r, t) + 2g|$(r, t)\'2p\r, t)ij{r, t) 



+ 



|<l>2(r, t)^t(r, t)^t(r, t) + |$* ^(r, t)^(r, t)^{r, t) 
d^r ip^r, t)4}\r, t)Vi(r, t)i>{r, t) , 



(18) 



where we have indicated with K2 the term quadratic in ip, ip"^, while K4 is the term quartic 
in these operators. By expanding K2 up to terms linear in the fluctuations (5$(r, t) we can 
rewrite it as the sum K2 = Klf^ + k!^^^ , where K^^ does not contain the fluctuations of the 
condensate 



d-'v 



2m 



+ Ve^t{r) - l^ + 2gno{r) ]'^{r,t) 



+ |no(r) {ij^{r,t)i,^{T,t)+ij{r,t)i>{T,t) 



(19) 



while is hnear in the fluctuations 5$ 



2g%(v) 5$(r, t) + 5^*{v, t) W{v, t)^/j{r, t) 



+ g<l>o(r) (5<^{v, t)i,\v, t)i,\T, t) + 5$*(r, t)V^(r, t)i,{v, t) jj . (20) 

In the quartic term we first use the mean field decomposition (all quantities depend on 
r and t) 

■ip'^tp^^^ = 4:htp^^ + fh^'^tp^ + m*'(p^ , (21) 

and then we expand the resulting expression up to linear terms in the fluctuations 5n and 
Srh thereby obtaining 



A5h{v,t)i)\v,t)i){v,t) 



+ ( (5m(r, t)V^^(r, t)ip\v, t) + 5m* {r, t)ip{v, t)ijj{r, t) 



(22) 



where we have made use of the Popov ansatz (|T2). In the above equation Kf"* is the zeroth 
order term, while K^l"* is linear in the fluctuations 5h and 5rh. In the equations of motion 
(P^) the commutators of a\aj and a^Uj with k'^'' yield, in the small-amplitude regime, the 
coupling to the fluctuations of the condensate, whereas the commutators with K^'^ give the 
coupling to the fluctuations of the normal and anomalous particle densities. If the density 
of the non-condensate particles is much smaller than the density tiq of the condensate, the 
coupling to the condensate is more important than the coupling to 5h and 5m and we can 
consequently neglect the contributions to the commutators arising from the term K^^^ . 

One can easily show that the operator i^^2°'* + -^4°^ is diagonal in the quasiparticle oper- 
ators aj, a^j if the functions Uj and Vj satisfy the coupled Bogoliubov equations 

Cvj{r) + gno{r)uj{r) = -ejVj{r) , (23) 
where we have introduced the hermitian operator 
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2m 



+ Kx-t(r) -/i + 2gn(r) 



As a consequence the relevant terms in the Hamiltonian become 

K2 + K^ = J2eja]aj+Ki^^ 



^2 5 



(24) 



(25) 



where the quasiparticle energies ej are obtained from the solutions of the Bogoliubov equa- 
tions 

The commutators in the equations of motion ([T7|) can be now calculated straightfor- 
wardly. To lowest order in the fluctuations one gets 
d 



^^QlU^) = (9 - + 2g(/° - /;) / <|.o(r) 



5$(r,t) + 5<l>*(r,t)jx 
u,(v)u*(v) + Vi{v)v*{j)^ + 5^Y,t)vi{r)u*{v) + 5^*{Y,t)Ui{Y)v*{v) 
for the time evolution of /jj, while the equation of motion for gij is given by 



(26) 



^^Ql9^J{t) = [e, + e,)g,,{t) + 2g(l + /« + /°) j %{r) 



<(rX(r) +<(r)M*(r^ 



5$(r, t)M*(r)M*(r) + 5$*(r, t)<(r)t;*(r) 



(27) 



In the equations (p6D-(p^) fj is the equilibrium density of quasiparticles 

/° = = (exp(/5e,) - 1)-^ (28) 

in terms of which the non-condensate particle density, at equilibrium, is written as 



E[(h-(r)r + |.;,(r)|V° + h-(r)p] . (29) 

j 

The last term in (^) accounts at T = for the quantum depletion of the condensate. The 
fluctuations 6n and 6m of the normal and anomalous particle densities can be straightfor- 
wardly expressed in terms of fijit) and Qijit) and the equation (|13D for the oscillations of 
the condensate takes the flnal form 



d 



2V72 



rv 

2m 



Kxi(r) + 2gn(r) 5$(r, t) + gno(r)5$*(r, t) 



+ g$o(r)E 



2 <(r)n,(r) + <(r)i;,(r) +<(r)«,(r) /,,(t) 



(30) 



+ 2t;,(r)M,(r) + w,(r)«,(r) b,,(t) + 2<(r)t;;(r) + v*{r)v*{r) )gUt) 
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The equations of motion (p6D-(p7D and (^) describe the small- amphtude coupled oscillations 
of the condensate and non-condensate particles in the coUisionless regime. They represent 
the main result of this section. 



III. DAMPING OF THE OSCILLATIONS 



Let us suppose that the condensate oscillates with frequency ui 



(31) 



Through the coupling term in the equations of motion (|26|)- (|27|) the fluctuations of the 
condensate act as a time-dependent external drive inducing oscillations in fij and Qij. The 
Fourier component of fij at the driving frequency u is given by 



huj + (et — e,) + iO 



5$1 [UiU* + ViV* + ViU* 



+ 5$2 (uiU* + ViV* + UiVj'j 



(32) 



where the frequency u has been chosen with an infinitesimally small component on the 
positive imaginary axis. Analogously, for the component of Qij oscillating at the frequency 
u one finds 



gij{uj) = 2g 



1 + /° + /° 

hu — (ej + ej) + iO 



d?Y <l>r 



5$! [u*v*+v*u*+u*u] 

+ 5$2 {ulv* + Vlu* + vlv*) 



(33) 



If one neglects in equation (pOD the coupling terms involving the fiuctuations of the non- 
condensate particles, the condensate components 5^\, 5$2 satisfy the unperturbed RPA 
equation |H] 



(34) 



/ 


c 


gno(r) ^ 


^(5$?(r) 






^5$?(r) 


\ 








= TiIjJq 




v 


gno(r) 




^<5$°(r) 


} 






) 



and the normalization condition 
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0|2 
2I 



(35) 



where uq is the unperturbed RPA eigenfrequency of the mode. We treat the couphng terms 
in equation (|30D as small perturbations and we write the solution of the normal mode as 



y5$2 J 



and the perturbed eigenfrequency as 



(36) 



hu = hujQ + 6E — i'j . 



(37) 



In eq. (0) 6E represents the shift in the real part of the frequency and 7 is the damping 
coefficient, while the eigenvector correction in eq. ( p6D is chosen to be orthogonal to the 
unperturbed eigenvector 



(38) 



Starting from the equations for the perturbed components and 6^2 [eq- 
complex conjugate] we multiply the first equation by 5$* and the latter by (5$2) take the 
difference of the two equations and finally integrate over space. By using (^) and (^) we 
get the following relation for the perturbed eigenfrequency 



I4-.|2 



+ 2g^E(i + /r + /;)' 



l-tJjj I 



,j ^hujQ - {ei + ej) + iO hujo + [ei + ej) + iO ^ 

where the matrix elements Aij, Bij and Bij are respectively given by 



(39) 



A 



B 



B; 



d^r $0 
d^r $0 
d^r $0 



5$?(nin* + ViV* + ViU*) + 5$^(uiU* + ViV* + Uiv* 
6^i{uiVj + ViUj + ViVj) + 5^^{uiVj + ViUj + UiUj) 



(40) 



The imaginary part of the right hand side of eq. (^) gives the damping coefficient 7. 
There are two distinct contributions to 7: one arises from the process of one quantum of 
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oscillation huo being absorbed by a thermal excitation with energy which is turned into 
another thermal excitation with energy = + hujQ. This mechanism is known as Landau 
(L) damping and is given by the imaginary part of the second term on the right hand side 
ofeq. (H) 

7L = 47r/ Y: \A,\' (/° - /°) Sihuo + e, - e,) . (41) 

ij 



The above result coincides with the finding of Ref. ||TJ] obtained by direct use of perturbation 
theory. A different mechanism of damping comes from the process where a quantum of 
oscillation huJo is absorbed and two excitations with energies + ej = hujQ are created. The 
damping associated with the decay of an elementary excitation into a pair of excitations, 
which has been studied by Beliaev in the case of uniform Bose superfiuids [^], is present also 



at T = 0, but is not active for the lowest energy modes in the case of a trapping potential 
because of the discretization of levels. The Beliaev [B) damping is given by the imaginary 
part of the first term in brackets on the right hand side of eq. (^) and reads 

7B = '2^9^ E \B^3? (l + + /°) ^(^^0 - Q - e,) . (42) 

The total damping coefficient is the sum of the two contributions 7 = 7l + 7_b, the Beliaev 
damping becomes dominant at low temperatures (fc^T <^ huo), while the Landau damping 
becomes dominant in the opposite regime of temperatures {ksT ^ huJo). 

In the case of magnetically trapped Bose gases result (^1]) and (^21) give the damping 
coefficient 7 both at low and high temperatures in the collisionless regime. This quantity 
has been measured at JILA over a wide range of temperatures for the m = and m = 2 
modes 0, revealing a very strong temperature dependence. To carry out the calculation 
at a given temperature T one must start from the unperturbed condensate eigenfrequency 
with the required symmetry, obtained from eq. (|3^), and the elementary excitation energies 
ej calculated from eqs. (pSf). One has then to calculate the matrix elements Aij and Bij 
(see eq. (^)) and then carry out the double summation in (^) and (^^. The explicit 
calculation of 7 in magnetically trapped gases will be the object of a future work. 
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By analyzing the real part of eq. (|39|) one can calculate the eigenfrequency shift 5E 
in the Popov approximation. However, for temperatures ksT < fi the effects arising from 
the equilibrium anomalous density m'^(r) are important for the calculation of the frequency 
shift and can not be neglected. In section 4-C we will present a calculation of the velocity 
of sound at T = for a uniform gas and it will be shown that the inclusion of m° is crucial 
to obtain the correct result. From this finding we conclude that a reliable calculation of the 
frequency shift 6E should be based on a more accurate dynamic theory, beyond the simple 
Popov ansatz (|Tll). 

IV. UNIFORM BOSE-CONDENSED GASES 



In this section we apply the results (|41| ) -(|42| ) to uniform Bose gases for which we reproduce 
well known results for the damping of phonons, both at finite and zero temperature, in the 
collisionless regime. 

For homogeneous systems the stationary condensate wavefunction is constant throughout 
space $0 = ^/no, while the condensate fluctuations and the excitations are described by plane 
wave functions 

7/ 7; 7/ 7; 

5$i(r) = -^e^^- 6M^) = -^e^""-^ «p(r) = -^e^P- Vp{r) = -^e^^^ , (43) 
where Up and Vp are real functions defined through the usual Bogoliubov relations 

^v = ^ + -P = 

and ep is the energy of the elementary excitations as obtained from the Bogoliubov equations 



2 \ 2 



1/2 



In these equations no = no(T) is the condensate density at the given temperature T obtained 
from the self-consistent calculation in thermodynamic equilibrium (see Ref. 0). 
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A. Thermal regime hujQ <ti ksT 



We consider the long-wavelength limit for the fluctuations of the condensate 

hujo = eq c:^cq , (46) 



where c = ygno/m is the velocity of sound calculated at temperature T. In this limit the u 
and V functions associated with the oscillations of the condensate can be expanded as 



As it has been shown in Ref. [14], using the above expansion for Ug and Vg one obtains the 
following low-g behaviour for the matrix elements App/ defined in (^) 

- (^) " + 4 + - ^-«|||) . (48) 

where 6ppr is the Kronecker symbol, 6 is the angle formed between the directions of p and q, 
and Vg = dep/ dp is the group velocity of the excitations. By introducing the above expansion 
into eq. ( PI) and carrying out the integration over 9 one finds the result (see Ref. [Q) 

l:,2L^L^„S'^F^r) , (49) 

where r = ksT /mc^ is a reduced temperature and F{t) is the dimensionless function 

FW=8Vj/rf.(e'-e-)-^(l-i--^)' . (50) 



where we have introduced the quantity u = Vl + ^t'^x'^. 

For temperatures r ^ 1 the function F takes the asymptotic limit F 37r^/^r/4 and 
the damping coefficient is given by 

7 SnkBTa . 

7g = ■ (^^^ 
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This regime of temperatures has been first investigated by Szepfalusy and Kondor As 
discussed in Ref. the numerical coefficient in (^) coincides with the one obtained in 
pO| , while it slightly differs from the one of p3|] . 

In the opposite regime of temperatures, r <^ 1, one finds the asymptotic limit F — *■ 
/5 and one recovers the well known result for the phonon damping |T^,|2^J2^,|T^ 



7 37r3 {kBTf 



8 mnTc'd' 



(52) 



In Fig. 1 the dimensionless function F{t) is plotted as a function of r together with its 
asymptotic behaviour both at small and large r's. One can see that F departs rather soon 
from the low-temperature r"^ behaviour, while it approaches the high-temperature linear law 
very slowly. 



B. Quantum regime %ujq ^ ksT 

At T = the damping of the long-wavelength fluctuations of the condensate with energy 
(^) is obtained through result with = 0. For uniform gases the matrix elements 
Bpp' entering eq. P^ ) read 



Bpp' — 5p' _p+q 



[uq{UpVp-q + VpUp-q + UpUp-g) + Vg{UpVp-g + VpUp-q + VpVp-q)^ . (53) 



'V 

In the Beliaev damping mechanism the momenta of the three excitations involved in the 
process are comparable, g ~ p ~ |p — q|- For p <^ mc one can use the following expansions 
for the excitation energies ep and the functions Up and Vp 

epC^ cp + 



„^ ^ H. 1 l^f H. i (JL\"" - i i^)"'' (54) 

^ \2p I 2\2mcJ 8\2mcJ 8 \2mcJ ^ ' 



mc 



1/2 1 / p X 1/2 I . p X 3/2 I . p x5/2 

^~~l2pi ^2VhncJ ~ 8VhncJ ~ 8\2mcJ 



If one substitutes the above expressions for all the m's and f 's in ( [53| ) and makes use of the 
condition for energy conservation 
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p-q|+p = g + ^^(9'+p' + |p-q|') , (55) 



after long but straightforward algebra one gets the result 



R -A 3 (gp|p-q|)V^ 



Once the low-energy behaviour of the matrix elements of B has been obtained, the calcula- 
tion of the damping coefficient follows directly 

1-1B = 7--^3 ' (57) 

and coincides with the well known result first obtained by Beliaev using diagrammatic 
techniques p^j26[| (for a review of the second-order Beliaev approximation at T = and its 



extension to finite temperature see Ref. pO|| ) 



C. Bulk compressibility 

In this section we calculate the isothermal inverse compressibility = n^[d^/dn)v,T- 
In a uniform system at T = it is fixed by the velocity of sound through the relation 

^T=^ = p{T] =pc'iT = 0), (58) 

where p = mn is the mass density. We will show that the Popov approximation gives an 
incorrect result for the inverse compressibility in the low-temperature regime and that one 
must go beyond this approximation in order to obtain the low-T behaviour of k^^. 

Without making the Popov assumption (pH]) and including the term proportional to 
m^{r), the stationary equation for the real wavefunction $o(r) becomes 



2m 



K.t(r) -fi + g(no(r) + 2n°(r) + m%r)) U'o(r) = , (59) 



instead of eq. ([T^). In terms of the quasiparticles operators defined in (|T^) the anomalous 
density at equilibrium is given by 
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m°(r) = Y: 2u,{v)v*iv)f^ + u,iv)v*iv) . (60) 

j 

For a uniform system the above equation (|59D fixes the value of the chemical potential 

= gn + g{h^ + fhP) (61) 

which differs from the value /xp = g{n + -nP) obtained in the Popov approximation At 
T = this approximation gives /xp(T = 0) = ATcfi^an/m{l + 8{a^ny/'^ /3^/7l), where the 
correction to gn comes from the quantum depletion of the condensate. The a^/^ order of the 
correction is correct but the numerical coefficient is wrong. To obtain the T = chemical 
potential to order a^^^ one must use eq. (0) and expand the coupling constant g to second 
order in the scattering length 

m \ m V p p / 
The need of the above renormalization of g for the calculation of the ground-state energy 



of a Bose gas, has been first pointed out by Lee, Huang and Yang and is discussed in 
many textbooks (see e.g. Ref. |2^). By substituting the renormalized value of g into the 
first term of eq. ( pT] ) one gets 

/m 4:iTh^an ( ^ 8 / m/2 47r/i^a I ( ni 1 \\ 

m \ 3V7r ) 

which coincides with the well-known result for the chemical potential of a dilute Bose gas at 
T = In the above equation we have used (^)-(|4^) for the product Mp(r)fp(r) with 
the Bogoliubov spectrum (^51) . Notice that in the integral over p the ultraviolet divergencies 
arising from the renormalization of g and from fhP cancel out. By differentiating (p3D with 
respect to n one obtains for the velocity of sound 



e(T = 0) = ,l^fl + 4(a3„)vA , (64) 



m? \ \/T[ 
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a result which has been first derived by Behaev [^. In a self-consistent dynamic theory 



result ( [6^ ) should also be obtained from the frequency shift of the long- wavelength excita- 
tions. However, it can not be found within the Popov approximation because the crucial 
ingredients, renormalization of g and equilibrium value of the anomalous density, are not 
accounted for in this approximation. 
At finite T one gets from eq. (|6l|) 

Anh^an ( 32 no 3 i/2\ 

/i T = 1 + —= — [a no ' 

m \ Sa/tt n J 



H — (a no) ' r / ax , , (65) 

where r = kBT/mc^ is the reduced temperature. Result (^) coincides with the finding of 
Ref. obtained using the second-order Beliaev approximation. In the low-temperature 



regime, r -C 1, one has 



Mr)^,(T^O) + ^(^ (66) 



for the chemical potential, and 

=pc (T = 0)---^ (67) 

for the inverse compressibility. The above results coincide with the ones obtained from the 
thermodynamic relation fi = [OFpfJ dN)yT-, where Fp^ is the free energy of a phonon gas. 
In the same temperature regime, fc^T -C mc^, the Popov approximation gives /ip(T) = 
/ip(T = 0) + m?c{kBTY / {\2nTt') and the inverse compressibility exhibits an unphysical 
dependence which is not consistent with the dependence obtained by differentiating the 
phonon free energy. 

In the low-temperature limit the frequency shift of coUisionless phonons is proportional 
to T^log(T) p^ , p5| , |3T| , whereas the Popov approximation again yields an unphysical 



dependence arising from the low-T expansion of no(T). The analogy with the result for the 
chemical potential suggests that also for this calculation the inclusion of fhP is crucial to 
obtain the correct result. 
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From the above results we conclude that a self-consistent dynamic theory, aiming to 
describe both the damping and the frequency shifts of the oscillations of the condensate, 
should go beyond the Popov ansatz (|11|). The ingredients that this theory should contain 
are the following: i) the equilibrium anomalous density m° has to be taken into account, ii) 
the renormalization ( p^ ) of the interaction coupling constant is needed in order to reproduce 
the energetics at T = 0, iii) the elementary excitation energies have to be gap less and 
must coincide with the Bogoliubov spectrum at low-temperatures. 
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FIGURES 

FIG. 1. Dimensionless function F as a function of r (solid line). The asymptotic behaviours 
for T <C 1 (dashed line), and for r 3> 1 (long-dashed line) are also reported. 
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